Statistics of opinion domains of the majority-vote model on a square lattice 
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The existence of juxtaposed regions of distinct cultures in spite of the fact that people's beliefs 
have a tendency to become more similar to each other's as the individuals interact repeatedly is a 
puzzling phenomenon in the social sciences. Here we study an extreme version of the frequency- 
dependent bias model of social influence in which an individual adopts the opinion shared by the 
majority of the members of its extended neighborhood, which includes the individual itself. This is 
a variant of the majority-vote model in which the individual retains its opinion in case there is a 
tie among the neighbors' opinions. We assume that the individuals are fixed in the sites of a square 
lattice of linear size L and that they interact with their nearest neighbors only. Within a mean-field 
framework, we derive the equations of motion for the density of individuals adopting a particular 
opinion in the single-site and pair approximations. Although the single-site approximation predicts a 
single opinion domain that takes over the entire lattice, the pair approximation yields a qualitatively 
correct picture with the coexistence of different opinion domains and a strong dependence on the 
initial conditions. Extensive Monte Carlo simulations indicate the existence of a rich distribution 
of opinion domains or clusters, the number of which grows with L 2 whereas the size of the largest 
cluster grows with InL 2 . The analysis of the sizes of the opinion domains shows that they obey a 
power-law distribution for not too large sizes but that they are exponentially distributed in the limit 
of very large clusters. In addition, similarly to other well-known social influence model - Axelrod's 
model - we found that these opinion domains are unstable to the effect of a thermal-like noise. 

PACS numbers: 89.75.Fb, 87.23.Ge, 05.50.+q 



I. INTRODUCTION 



Despite the tendency of social actors to become more 
similar to each other through their local interactions 
(HQ, global polarization, i.e., the existence of stable mul- 
ticultural regimes, is a ubiquitous phenomenon in human 
society 0, 0] • This issue has been addressed by somewhat 
idealized agent-based simulations of human behavior pro- 
ducing highly nontrivial insights on the nature of the col- 
lective behavior that results from the homogenizing local 
interactions (see Q for a recent review). 

In this line, a particularly successful model is Axel- 
rod's model for the dissemination of culture or social in- 
fluence @. In Axelrod's model, the agents are placed 
at the sites of a square lattice of size L x L and can 
interact with their nearest neighbors only. The culture 
of an agent is represented by a string of F cultural fea- 
tures, where each feature can adopt a certain number q of 
distinct traits. The interaction between any two neigh- 
boring agents takes place with probability proportional 
to the number of traits they have in common. Although 
the result of such interaction is the increase of the simi- 
larity between the two agents, as one of them modifies a 
previously distinct trait to match that of its partner, the 
model exhibits global polarization Q. The key ingredi- 
ent for the existence of stable globally polarized states in 
Axelrod's model is the rule that prohibits the interaction 
between agents which do not have any cultural trait in 
common. Relaxation of this rule so as to permit interac- 
tions regardless of the similarity between agents leads to 
one of the q F distinct absorbing homogeneous configura- 
tions 7]. In addition, introduction of external noise so 



that traits can change at random with some small prob- 
ability Q, as well as the increase of the connectivity of 
the agents (9l4TTI|. destabilizes the polarized state also. 

Although Axelrod's model enjoyed great popularity 
among the statistical physics community due mainly to 
the existence of a non-equilibrium phase transition (l2l | 
that separates the globally homogeneous from the glob- 
ally polarized regimes [lj, [wj, the vulnerability of the 
polarized absorbing configurations was considered a ma- 
jor drawback to explaining robust collective social be- 



havior. In this vein, Parisi et al. [15| have proposed a 
lattice version of the classic frequency bias mechanism 
for cultural or opinion change 0,13]) which assumes that 
the number of people holding an opinion is the key fac- 
tor for an agent to adopt that opinion, i.e., people have 
a tendency to espouse cultural traits that are more com- 
mon in their social environment. Actually, almost any 
model of cultural transmission admits that the probabil- 
ity that an individual acquires a cultural variant depends 
on the frequency of that variant in the population. The 
frequency-dependent bias mechanism requires, however, 
that the individual be disproportionately likely to acquire 
the more common variant [3|. 

More to the point, in the model of Parisi et al. the cul- 
ture of an agent is specified by a binary string of length 
F (so q = 2) and each bit of that string takes the value 
which is more common among its neighbors [l5| . Since 
the model can be immediately recognized as F indepen- 
dent majority- vote models [l|| we found the claim 
by those authors that such model exhibits a polarized 
regime most intriguing. In order to check whether the 
multicultural absorbing configurations reported by Parisi 
et al. were not artifacts of the small lattice size used in 
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their study (L = 20), in this paper we present a detailed 
analysis of the effects of the finite size of the lattice by 
considering square lattices of linear size up to L = 4000. 
In addition, we set F = 1 to stress the identification with 
the well-known majority vote model of statistical physics. 
Hence the agents or sites of the lattices are modeled by 
Ising spins. 



Our findings indicate that the polarized regime is in- 
deed stable in the thermodynamic limit L —> oo and that 
the element accountable for this stability is the proce- 
dure used to calculate the majority, which includes the 
site to be updated (target site) in addition to its near- 
est neighbors. For sites in the bulk of the lattice, this 
procedure essentially specifies the criterion of update in 
case of a tie, i.e., in case there is no majority among the 
neighbors. In this case, the majority-vote variant used 
by Parisi et al. leaves the state of the agent untouched, 
whereas the most common variant used in the literature 
sets that state at random with probability 1/2 [TH fl9|. 
It should be noted, however, that the criterion used by 
Parisi et al. is actually the original definition of the 
majority- vote model as introduced in Refs. [l6L fl7|. It is 
surprising that such (apparently) minor implementation 
detail produces nontrivial consequences in the thermody- 
namic limit, which can actually be confirmed analytically 
using a mean-field approach that takes into account the 
correlation between nearest neighbors - the pair approxi- 
mation [12] . However, the addition of a thermal-like noise 
such that the individuals may take a stance opposed to 
the majority opinion with some small probability desta- 
bilizes the globally polarized regime leading to one of 
the two low-temperature homogeneous steady states of 
the Ising model [20j . In that sense we disagree with the 
claim by Parisi et al. that their variant of the majority- 
vote model is more robust than Axelrod's model to the 
effect of this type of noise. 



The rest of this paper is organized as follows. In 
Sect. HJ we describe the variant of the majority-vote 
model proposed by Parisi et al., which henceforth we 
refer to as the extended majority- vote model since its 
key ingredient is the stretching of the neighborhood to 
include the target site [H)]. In Sect. IIIII we present 
the mean-field approach to this model using the single- 
site and pair approximations [H, 0|. In Sect. IIVI we 
use extensive Monte Carlo simulations to investigate sev- 
eral properties of the absorbing configurations such as 
the average number of opinion domains or clusters, the 
size of the largest cluster and the distribution of cluster 
sizes, giving emphasis to their dependences on the lattice 
size. In that section also we discuss the vulnerability of 
the heterogeneous opinion regime against a thermal-like 
noise, which allows the individuals to disagree with the 
majority opinion. Finally, in Sect. [V] we present our 
concluding remarks. 



II. MODEL 

The agents are fixed at the sites of a square lattice of 
size L x L with open boundary conditions (i.e., agents 
at the corners interact with two neighbors, agents at the 
sides with three, and agents in the bulk with four near- 
est neighbors). The initial configuration is random: the 
opinion of each agent is set by a random digit or 1 
with equal probability. At each time we pick a target 
agent at random and then verify which is the more fre- 
quent opinion (1 or 0) among its extended neighborhood, 
which includes the target agent itself. The opinion of the 
target agent is then changed to match the corresponding 
majority value. We note that there are no ties for agents 
in the bulk or at the corners of the square lattice since in 
these cases the extended neighborhood comprises 5 and 
3 sites, respectively. However, agents at the sides of the 
lattice have an extended neighborhood of 4 sites (i.e., 3 
neighbors plus the target agent) and so in case of a tie, 
the opinion of the target agent remains unchanged. Of 
course, in the limit of large lattices the contribution of 
these boundary sites will be negligible. This procedure 
is repeated until the system is frozen in an absorbing 
configuration. 

Although the majority-vote rule or, more generally, 
the frequency bias mechanism for opinion change Q is 
a homogenizing assumption by which the agents become 
more similar to each other, the above-described model 
does seem to exhibit global polarization, i.e., heteroge- 
neous absorbing configurations (l5j . Since the study of 
Ref. [HI was based on a small lattice of linear size L = 20 
a more careful analysis is necessary to confirm whether 
this conclusion holds in the thermodynamic limit as well. 
It is interesting to note that for the more popular vari- 
ant of the majority-vote model, in which the state of 
the target site is not included in the majority reckoning, 
and ties are decided by choosing the opinion of the target 
agent at random with probability 1/2, the only absorbing 
states in the thermodynamic limit are the two homoge- 
neous configurations. For finite lattices, however, we find 
heterogeneous absorbing configurations characterized by 
stripes that sweep the entire lattice. In Sect. IIVI we 
carry out a finite size scaling analysis of the extended 
majority-vote rule aiming at understanding the nature 
of the opinion domains that fragment the absorbing con- 
figurations. 



III. MEAN-FIELD ANALYSIS 

In this section we offer an analytical approximation to 
the extended majority-vote model introduced by Parisi 
et al 15[ . The state of the agent at site i of the square 
lattice is represented by the binary variable r\i = 0, 1 
and so the configuration of the entire lattice comprising 
N = L 2 sites is denoted by r\ = (771,772,...,^). The 
master equation that governs the time evolution of the 
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probability distribution P (r],t) is given by 
jP („, t) = [Wi (ff) P (fj\t) - Wi (v) p (v, *)] (1) 

i 

where rf = . . . , 1 — Vh . . . , Vn) and Wi (77) is the 
transition rate between configurations r\ and rf [l8|, [l!| . 
For the extended majority-vote model we have 







^ Vi+s +rji-3 
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(2) 



where the notation 5Z 5 (. . .) stands for the sum over the 
4 nearest neighbors of site i and O (x) = 1 if x > and 
otherwise. 

We are interested in determining the fraction of agents 
holding opinion 1 or, equivalently, the fraction of sites 
in state 1, which we denote by p. Since all sites are 
equivalent we have p = X^ 7 /*/-^ = (Vi) an d this mean 
value is given by the equation 

±( r)i ) = {(l-2r ]i )W i (r ] )) (3) 

with ((. . .)) = J2 V (• • •) P (Vt *) as usual. To carry out this 
average we need to make approximations. In what fol- 
lows we study in detail two such approximation schemes, 
namely, the single-site approximation and the pair ap- 
proximation. 



A. The single-site approximation 

This the simplest mean-field scheme which assumes 
that the sites are independent random variables so that 
Eq. ([3]) becomes 



dp 
~dt 



.£B*(p)|e[n-2]-l| 
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•(l-p)£iJ«(p)e[n-3] (4) 
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where B™ is the Binomial distribution 
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Carrying out the sums explicitly and rearranging the 
terms yield 

g = -p(l-p)(2p-l)(3p 2 -3 P -l). (6) 

This equation, which is invariant to the change p <H> 1 — p, 
has three fixed points, namely, p* — 0, p* — 1 and 
p* = 1/2. The first two fixed points are stable and the 
third one is unstable. This means that in the single- 
site approximation the only stable configurations are 
the homogeneous ones. Finally, we note that p con- 
tains the same information as the single-site probabil- 
ity distribution pi (rji). In fact, p\ (r/i = 1) = p and 
Pi (m = 0) = 1 - p- 



B. The pair approximation 

In this scheme we assume that nearest neighbors sites 
are statistically dependent random variables so that, in 
addition to the single-site probability distribution pi (rji), 
we need to compute the pair probability distribution 
Vi ("Hii Vi+s) as well. Of particular importance for the sub- 
sequent calculations is the conditional probability distri- 
bution pi\i (r/i + s | f]i) which is given simply by the ratio 
between the pair and the single-site probability distribu- 
tions. 

For the sake of concreteness, let us denote the states 
of the 4 neighbors of site i by 771, 772, r/3 an( i Vi- These are 
independent random variables in the pair approximation 
since they are not nearest neighbors in the square lattice, 
so the sum n = rji + V2 + f]3 + ^4 that appears in the 
argument of the Theta functions is a sum of independent 
variables. With these notations we can rewrite Eq. §3§ 
as 

^ = ~P P4|i07i>»72>»73>»74 I l)|0(n-2)-l| 

rjl,...,?74 

+ (!-/>) Y PMi('ni,V2,V3,m I 0)8(n-3) 



rjl,...,J)4 



(7) 



where we have carried out the sum over rji = 0, 1 ex- 
plicitly. The assumption that 771 , . . . , 774 are statistically 
independent random variables allows to write 

PA\l(Vl,V2,V3,V4 I Vi) =Pl\l(Vl I Vi) x ■■■ ><Pl\l(V4 I Vi) 

(8) 

and finally obtain 

^ = -p£#[pi|i(i|i)]ie[»-2]-i| 

71 = 

4 

+ (l-p)^B^[p 1 | 1 (l|0)]e[«-3] (9) 

n=0 

which is identical to Eq. Q except for the arguments 
of the binomial distributions. Using the notation <p = 
P2 (1, 1) we write 



Pi|i(l|l) = 
(1 I 0) = 



(10) 

(11) 



so that Eq. (O involves only two unknowns, p and <f>. 
Carrying out the summations explicitly so as to eliminate 
the Theta functions yields 
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FIG. 1: Time evolution of the fraction of sites in state 1, p, in 
the pair approximation obtained by solving Eqs. (|12[) and (|15[) 
using Euler's method with step-size 0.01 for different initial 
conditions (bottom to top) pa = 0.1, . . . , 0.9. 



We note that this equation reduces to Eq. © in the case 
the neighboring sites are assumed independent, i.e., <f> — 
P 2 . 



Q. 




FIG. 2: The fraction of sites in state 1 at equilibrium, p* 
(represented by o), and the probability that two neighbors are 
in state 1 at equilibrium, (j>* (represented by A) as functions 
of the initial fraction of sites in state 1, po- The solid line is 
the result of the pair approximation for which cf>* —tp*. The 
initial condition is p = po and <j) — p§. The symbols show the 
results of the Monte Carlo simulations for a square lattice of 
linear size L = 200. 

Next, our task is to determine an equation for <fi = 
{ViVj) where j labels one of the four nearest neighboring 
sites of i. We have 

j t {rnnj) = 2 ( Vj (1 - 2th) Wi (v)) ■ (13) 

Carrying out the average over rji and r\j (say j = 1) 



explicitly yields 

jj. 3 

' ' = -2 PPlll B l [Pi\i (1 I 1)] I© [m - 1] - 1| 



dt 



m=0 



-2 (1 - p) Plll (1 | 0) ]T Bf n [ Pl|1 (1 | 0)] 9 [m - 2] 

(14) 



m=0 



which can be written more compactly as 



ft = u 

+2 (p - <j>) 



B: 



p- 4> 
l-p 



p~<t> 
i-p 



(15) 

Equations (|12l) and (IT5|) determine completely the time 
evolution of p and <fi and so they are our final equations 
for the pair approximation of the extended majority- vote 
model. Figure [TJ which shows the time evolution of the 
density of sites in state 1, confirms that Eq. (fT5)) is indeed 
invariant to the change <H> 1. Most surprisingly, this 
figure uncovers an unexpected dependence on the choice 
of the (random) initial condition po an d (f>o = Po' indeed 
in the range po £ (p m , 1 — p m ), where p m ~ 0.25, the 
equilibrium solution p* is a smooth function of po, i.e., 
there is a continuum of fixed points. Accordingly, Fig. [2] 
shows the fixed point solutions of Eqs. (fl~2"|) and (fT5")) for 
the usual situation in which the states of the sites of the 
initial configuration are set 1 with probability po and 
with probability 1 — po- For this setting we have (po = Pq- 

The mathematical explanation for the odd dynamical 
behavior exhibited in Figs. [1] and [5] is that the dynamic 
evolution is such that p — > for t — > oo, a condition that 
solves the two equations (TT2l and (jl"5j) at equilibrium (i.e., 
dp/dt = d<fr/dt = 0) simultaneously and so leaves one of 
the unknowns free to take any arbitrary value set by the 
dynamics or the initial conditions. The physical reason 
is that there are in fact many absorbing configurations 
very close to the random initial configurations and a few 
lattice updates (typically 10) is sufficient to freeze the 
dynamics (see Fig. [3]). In fact, the comparison with the 
results of the Monte Carlo simulations exhibited in Fig. 
[5] shows a good qualitative agreement between the pair 
approximation and the simulation results. 

Finally, we note that the pair approximation for the 
variant of the majority-vote model in which ties are de- 
cided by flipping the state of the target site with probabil- 
ity 1/2 yields the anticipated result that the only stable 
fixed points are p = <j) — and p = cf> = 1, whose basins 
of attraction depend on whether po is less or greater than 
1/2. 



IV. MONTE CARLO SIMULATIONS 

In this section we first demonstrate that the heteroge- 
neous absorbing configurations of the extended majority- 
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FIG. 3: A typical absorbing configuration of the extended 
majority-vote model for a lattice of linear size L = 300. 
Agents with opinion 1 are painted black and agents with opin- 
ion are painted white. There are a total of M = 779 clusters 
and the largest one comprises S m ~ 4581 agents. 




L 2 



FIG. 4: Logarithmic plot of the average number of clusters 
(M) as function of the the lattice area L 2 . The solid straight 
line is the fitting (M) = 0.00742L 2 for large L. Each symbol 
represents the average over 10 4 samples so the sizes of the 
error bars are smaller than the symbol sizes. 



A. Statistics of the number of clusters 



vote model do persist in the thermodynamic limit and 
then we proceed with the characterization of the opinion 
domains (clusters) that fragment those configurations. 
Figure [3] illustrates one such a configuration for L = 300. 
We recall that a cluster is simply a connected, bounded 
lattice region wherein the agents share the same opin- 
ion. We focus on a few relevant statistical quantities to 
characterize the cluster organization, namely, the average 
number of clusters (M), the average size of the largest 
cluster (S m ), and the distribution of cluster sizes Ps- We 
must evaluate these quantities for different lattice sizes 
(typically we use 10 4 samples for each L) and then take 
an appropriate extrapolation procedure to infinite lat- 
tices [L — > oo). 

To simulate efficiently the extended majority-vote 
model for large lattices we first make a list of the active 
agents. An active agent is an agent whose opinion dif- 
fers from the most frequent opinion among its extended 
neighborhood. Clearly, only active agents can change 
their opinions and so it is more efficient to select the tar- 
get agent randomly from the list of active agents rather 
than from the entire lattice. In the case that the opinion 
of the target agent is modified by the majority-vote rule, 
we need to re-examine the active/inactive status of the 
target agent as well as of all its neighbors so as to update 
the list of active agents. The dynamics is frozen when 
the list of active agents is empty. The implementation of 
a similar procedure allowed the simulation of Axelrod's 
model for very large lattices and the clarification of sev- 
eral issues regarding the stability of the heterogeneous 
absorbing configurations of that model [3, [2l| . 



We begin our analysis by presenting the dependence of 
the average number of clusters (M) on the lattice area L 2 
in Fig. [H The increase of (M) with increasing L 2 is the 
evidence that confirms that the extended majority-vote 
model exhibits heterogeneous absorbing configurations in 
the thermodynamic limit, L — y oo. More to the point, 
we find 

lim = 0.00742 ± 10~ 5 . (16) 

For the purpose of comparison, for frozen random config- 
urations in which the L 2 sites are set to 1 or with the 
same probability this ratio yields 0.1342 ± 0.0001 (and 
so the average number of clusters scales linearly with the 
lattice area too). We note that the reciprocal of the ra- 
tio given in Eq. (TTBl is the average size of the clusters: 
(S) sa 134.8 for the extended majority-vote model and 
(S) « 7.45 for random configurations. 



B. Statistics of the largest cluster 

Although the mean size of a typical cluster is finite, 
there may be clusters whose sizes can become arbitrarily 
large as the lattice size increases. To investigate this 
point we show in Fig. [5] the average size of the largest 
cluster (S m ) as function of the lattice area. As in the case 
of (M), the asymptotic scaling is revealed only for large 
lattices (L > 500) and it indicates that (S m ) increases 
with InL 2 for large L. 

In the thermodynamic limit the relevant quantity is 
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FIG. 5: The average size of the largest cluster (S m ) as func- 
tion of the the lattice area L 2 . The solid line is the fitting 
(S m ) = 1933 In L 2 — 36363. Note the logarithmic scale in the 
x-axis. 



then 



lim = 1933 ± 10. 

L-i-oo In L 1 



(17) 



We found that (S m ) scales with \nL 2 for random con- 
figurations also, though the ratio liniL_ ) . 00 (5 m }/lnL 2 = 
53.0 ± 0.2 is considerably smaller than the ratio given in 
Eq. (HID. 

It is interesting to note that the standard deviation 

9 -|l/2 

(S 2 ) - (S) tends to the constant value 695.5 ± 0.5 

in the thermodynamic limit. This amounts to a large but 
finite variance and so in order to satisfy the Tchebycheff 
inequality the set of clusters whose size grows like lnL 2 
must have measure zero in that limit. Since the prob- 
ability that a randomly chosen site belongs to one such 
a cluster vanishes like (InL) /L 2 , the measure of the set 
of diverging clusters will be vanishingly small if its car- 
dinality is finite, i.e., if only a finite number of clusters 
have sizes scaling with In L 2 . 



C. Distribution of cluster sizes 

Up to now we found no qualitative differences between 
the statistical properties of the clusters produced by the 
extended majority- vote rule or by assigning randomly the 
digits 1 and to the sites of the square lattice. In fact, for 
both types of configurations the quantities (M) and (S m ) 
exhibit the same scaling behavior with the lattice area L 2 . 
However, a more detailed view of the cluster organization 
is given by the distribution Ps of cluster sizes S, which 
is shown in Fig. [5] for different lattice sizes. The data 
seems to be remarkably well fitted by the power-law dis- 
tribution Ps ~ S"" 1 ' 5 for over more than three decades. 
In addition, the figure indicates that the region where 



FIG. 6: Logarithmic plot of the distribution of cluster sizes 
S for (solid lines from left to right at Ps = 10 -6 ) L = 
50, 100, 200 and 1000. The dashed straight line is the fitting 
P a — 0.79S 1- ' 5 . These distributions are obtained by sampling 
10 7 absorbing configurations of the extended majority-vote 
model. 



the power-law fitting holds increases with increasing L. 
This power-law distribution, however, is not compatible 
with our previous findings of finite values for the mean 
and variance of Ps, since the distribution Ps ~ S* -1,5 has 
infinite mean and variance. 

To answer this conundrum, in Fig. [7] we plot the same 
data using a semi-logarithmic scale instead of the log- 
log scale of Fig. El The results show that for large lattice 
sizes the distribution Ps is in fact exponential, the power- 
law behavior being valid only in an intermediate range of 
cluster sizes. The results for the larger lattice L = 1000 
(though not the results for the smaller lattices L — 50 
and L = 100) show that the exponential distribution is 
not a mere cut-off inherent to the finitude of the lattices 
used in the simulations, as manifested by the straight-line 
appearance of Ps revealed by the semi-log graph in the 
regime of large clusters. For the sake of comparison, in 
Fig. [5] we show the distribution of cluster sizes for frozen 
random configurations which can be well-described by 
the exponential distribution P s w 0.00085 exp (-0.025) 
for large S. The puzzling power-law behavior found for 
intermediate values of S in the extended majority-vote 
model is missing in this case. 

Finally, it is interesting to note that, for any absorbing 
configuration of the extended majority-vote model, clus- 
ters of size 1 are prohibited, clusters of sizes 2 and 3 are 
allowed at the sides of the square lattice only, and the 
only cluster of size 4 allowed in the bulk of the lattice is 
that in which the four sites surrounding a unit cell are 
in the same state. In fact, this elementary square cluster 
is the only structure capable to propagate through the 
lattice, i.e., whenever a stable cluster is formed one such 
elementary square is formed. 
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FIG. 7: The same data shown Fig. \6\ plotted using a semi- 
logarithmic scale. 
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FIG. 8: Semi-logarithmic plot of the distribution of cluster 
sizes S for random configurations with (solid lines from bot- 
tom to top) L = 100 and 300. These distributions are ob- 
tained through the sampling of 10 8 random configurations. 
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FIG. 9: Fraction of sites in state 1 as function of the number 
of Monte Carlos steps r for the thermal noise intensity (solid 
lines from left to right) p — 10~ 3 , 10 -4 and 10 -5 and for a 
lattice of size L = 50. The dashed line is the evolution for 
the noiseless case which rapidly gets stuck in an absorbing 
configuration. The initial configuration is the same for the 
four trajectories. 

L = 50 we present in Fig. the time evolution of the 
fraction of sites in state 1, which we denote by p as in 
Sect. IIII1 for different values of the noise parameter p 
but for the same initial configuration. This figure, which 
exhibits the evolution of a single sample, shows the vul- 
nerability of the heterogeneous absorbing configurations 
against a vanishingly amount of thermal noise if one gives 
enough time for thermalization. For the small values of p 
considered, the steady-state is close to one of the two ho- 
mogeneous configurations, with all individuals exhibiting 
opinion or opinion 1. The specific outcome is a stochas- 
tic process which depends on the choice of the initial con- 
figuration as well as on the sequence of random numbers 
used in the update rule. 

V. CONCLUSION 



D. Vulnerability to noise 

A word is in order about the robustness of the hetero- 
geneous absorbing configurations to the effect of noise, 
which allows the agents to oppose the opinion of the ma- 
jority with some small probability. More pointedly, after 
applying the update rule for each agent, which may or 
may not result in the change of its opinion, we flip its 
opinion with probability p s [0,1/2] 19]. Since in this 
case there are no absorbing configurations the notion of 
active sites is useless and so we implement the dynamics 
by picking lattice sites at random with the same proba- 
bility. Then the unit of time (a Monte Carlo step) cor- 
responds to the update of L 2 randomly chosen sites. For 



This paper is motivated by the claim that the extended 
majority-vote model, for which the target site counts in 
the reckoning of the majority, exhibits a nontrivial multi- 
cultural regime, i.e., heterogeneous absorbing configura- 
tions fl5| . From the perspective of the statistical physics, 
we find this claim most intriguing, as we expect the 
model to behave similarly to the majority-vote models 
commonly considered in the physics literature, which are 
basically Ising models with zero-temperature Glauber ki- 
netics [20] , and so exhibit only two homogeneous absorb- 
ing configurations in the thermodynamic limit. Our first 
impression was that the heterogeneous absorbing config- 
urations were artifacts of the small lattice size (L = 20) 
used or of some 'non-physical' elements of the original 
model such as the nearest and next-nearest neighbors in- 
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teractions (Moore neighborhood) and the parallel update 
of the sites [l5[ . 

Accordingly, we have modified the original model pro- 
posed by Parisi et al. by introducing the usual near- 
est neighbors interactions (von Neumann neighborhood) 
and the random sequential update of sites. A careful fi- 
nite size scaling analysis of a few relevant measures that 
characterize the absorbing configurations shows that the 
heterogeneous absorbing configurations not only persist 
in the thermodynamic limit but produce a tricky distri- 
bution of cluster sizes Ps that decays as a power law 
Ps ~ S' 1 - 5 for clusters of intermediate size S and as 
an exponential for very large clusters sizes (see Figs. [5] 
and O . Essentially, the reason for this somewhat unex- 
pected outcome is the criterion of update in case of a tie 
for a site in the bulk of the lattice: the site remains un- 
changed rather than flipping to another state with proba- 
bility! /2 as usually done in statistical mechanics models 
[ll . [l9| . What is remarkable is that the existence of these 
heterogeneous absorbing configurations for the extended 
majority-vote model, as well as their absence in the more 
usual variants of that model, can actually be predicted 
analytically using the pair approximation: the multiple- 
clusters regime is signaled by the presence of an infinity 
of attractive fixed points of the mean-field equations. 

We found that the statistical properties (e.g., average 
number of clusters, mean size of the largest cluster and 
the asymptotic distribution of cluster sizes) of the clus- 
ters that break up the absorbing configurations (see Fig. 
[3]) are qualitatively identical to those of random config- 
urations. This could also be inferred by the very short 
convergence times (see dashed line in Fig. O which in- 
dicate that the absorbing configurations are very close 
to the initial random configurations, since a few lattice 
updates suffice to reach them. We are not aware of any 
similar study of the statistical properties of the hetero- 
geneous configurations of Axelrod's model, so we cannot 
compare the cluster organization of these two models. 



Our findings, which corroborate the results of Parisi 
et al. [H| , shows that the extended majority- vote model 
does exhibit a 'multicultural' regime, despite the fact that 
the update rule biases the agents to become more sim- 
ilar to each other. A similar conclusion holds for Ax- 
elrod's model as well @, [l3|, [l4|, except that in Axel- 
rod's model the similarity is a prerequisite for interac- 
tions - the 'birds of a feather flock together' hypothe- 
sis which states that individuals who are similar to each 
other are more likely to interact and then become even 
more similar (A similar assumption has been used 
to model the interspecies interactions in spin-glass like 
model ecosystem |22|.) The majority- vote model is con- 
siderably simpler and converges to the absorbing con- 
figurations much faster than Axelrod's. However, the 
(short) inventory of advantages stops here: in disagree- 
ment with the claim of Parisi et al. [15] we found that 
the absorbing configurations of the extended majority- 
vote model are vulnerable to the noisy effect of flipping 
the opinions of the agents with some small probability. 
It is well-known that this type of noise destabilizes the 
heterogeneous configurations of Axelrod's model too Q. 
Of course, the extended majority-vote model lacks the 
main appealing feature of Axelrod's model, namely, the 
existence of a non-equilibrium phase transition that sep- 
arates the homogeneous and the polarized regimes in the 
thermodynamic limit. In that sense it would be interest- 
ing to find out how a similar transition could be induced 
in the zero-temperature majority-vote model. 
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